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Abstract 

Two types of finite series of products of harmonic numbers involving nonneg- 
ative integer powers are evaluated, also yielding two other important harmonic 
number identities. The recursion formulas for these sums are derived, which are 
easily translated into a computer program. 
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1 Definitions and Basic Identities 

The generaUzed harmonic numbers used in this paper are: 

Ht' = E ^ (1-1) 

k=l 

from which follows that iJQ™' ~ 0. The traditional harmonic numbers are: 

Hn - H^'^ (1.2) 

A well known identity is [2 [31 El E] : 

n 

j:lH,^liH'^ + Hi'^) (1.3) 
fc=i 

and [an]: 

n 
fc=0 

Let be the Bernoulli number with = 1/2 instead of Bi = —1/2: 

S+ = (-1)"B„ =B„ + ^„,i (1.5) 
Then for nonnegative integer p [4] : 



and for nonnegative integer p [1]: 

n 



k=0 



^ fe=i 



p+1 
k 



(1.7) 



When a < b are two integers and {xk} and {yk} are two sequences of complex 
numbers, and {s^} the sequence of complex numbers defined by: 



Sk 



then there is the following summation by parts formula [1]: 



&-1 



b-l 



^ XkVk = Sb-iVb - ^ sk{yk+i ~ Vk) 



(1.8) 



(1.9) 



2 First Type of Harmonic Number Identities 

Theorem 2.1. For nonnegative integer n: 

n 

^ HkHr^-k = (n + l)[(-ff„+i - 1)' - H^^l^ + 1] (2.1) 



k=0 



Proof. With generating functions: let [z"]/(z) be the coefficient of z" in the power 
series expansion of /(z). Then [2]: 



1 — z 1 — z 



(2.2) 



The generating function of a convolution of two series is the product of their generating 
functions: 



$:M„-. = n(^^(hr(^)) 



fc=o 



These generating functions are known [2]: 

. 1 



(1 - zf 



[z"](ln( — )) 



1 - z' 



(2.3) 



(2.4) 



(2.5) 



where [^] is the Stirling number of the first kind, which in this case evaluates to [5]: 

« (2.6) 

(71 - l)!ff„_i if n > 



n 




2 





2 



The product of these two generating functions results in another convolution; 



11" 1 



(2.7) 



k = l 



which using (1.71 with p — and m — 1 and (1.4) results in: 

n-l 



fc=0 



fe=0 



= {n+ 1)[H^ - H^f>] - 2[nHn - n] 
With iJ,!™^ = h':^\ - l/(n + 1)'" this resuhs in the theorem. 



(2.8) 



□ 



Using different methods this theorem was already proved in literature [6l [7] . Now 
evaluating the sum directly yields: 

n n— 1 k n~k 



n— 1 ^ n — 1 n — fc ^ 

n— 1 ^ n— z ^ n— j 

i— 1 j — 1 A;— i 

n— 1 ^ n— -i 

i=i j=i 

n — l ^ n— 1 

= (n + l)^^i/„_, -^i/„_ 



(2.9) 



1=1 

ri-l 



1=1 
ri-1 



E 

i=l 



= (ji + 1) XI - E " "^"-1 + (" " 1) 



The second sum is evaluated using (1.7 1 with p ~Q and m = 1, and substituting the 



right side of (2.8 1 results in the following identity: 

" 1 



(2.10) 



fe=i 
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Evaluating the sum for n + 1 and using Hn^^+i = Hn-k + — A; + 1) yields: 

n+l n 



k=0 



k=0 



n n ^ 

fc=0 fe=0 



Hk 



(2.11) 



k=0 



k=Q 



Substituting (2.1) for n + 1 and n results in the following identity: 



(2.12) 



fc=0 



Using summation by parts (1.9) with Xk = Hn^k and yk — k^Hk, the following 
recursion formula results: 

Y k^HkH^-k = -^{{n + min + l)Hl ^{n- l)i/„ - ^] 
^-^ P + 1 n + 1 



k=Q 



n 

(n + l)i/„ ^[((fc + I)'' - k^Hk + {k + ly-'] 



k=0 



^[(n - k){{k + ir - F) + pkP]HkHr,- 

k=0 
n 

Y + lT-kP + {n-k + l)P-^]Hk 

k=0 
n 

j2k{k + ir-'} 



(2.13) 



fe=0 



For p = this formula yields (2.1) using (2.12). Letting R{n,0) be (2.1) and using 
F{n,p,m) from (1.7) evaluation of this formula results in the following recursion 
formula for integer p > 0: 



Yk''HkHn-k = R{n,p) 
1 



fc=0 



p + I 



fe=0 



p - k + 1 



)R{n,k) 



(2.14) 



+ (1 + (-l)^^-^(n + l)P-''-^)F{n, fc + 1, 1) 
P 
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For a list of the resulting expressions up to p = 5 see section Q below. 

3 Second Type of Harmonic Number Identities 

Theorem 3.1. For nonnegative integer n: 

n 

^ Hi = (n + - (2n + 3)il„+i + 2(n + 1) 



(3.1) 



k=0 



Proof. 



k=0 



n n k ^ 

71 ^ n 

i—l k—i 
n ^ n i—l 

= E7(E^^-E^^) 

1=1 k=l k = l 

n 

= H„[{n + l)ff„+i - (n + 1)] - V ^{iH, - ^) 

Z 1 



(3.2) 



i=l 



Hn[{n + l)Hn+i -{n + 1)] - (n + l)iJ„+i + 2n + 1 
{n + l)i/2^i - (2n + 3)i/„+i + 2(n + 1) 



□ 



Using summation by parts (1.9) with Xk = Hk and yu = k^H^, the following 
recursion formula results: 



1 



-{{n+l)P[{n+l)H^-{n-l)Hr, 



k=0 



n 

Y^[{k + i){{k + iY-kn-pkP]Hl 



k=0 



(3.3) 



j2[ik + ir-kiik+ir-kn]Hk 

n 

Y.k{k + \Y-^} 



k=0 



For p = this formula yields (3.1 ). Letting S{n, 0) be (3.1 1 and using F{n,p, m) from 
(1.7 1 evaluation of this formula results in the following recursion formula for integer 



5 



p > 0: 



fe=0 



-{(n+ l)[(n + - i/„+i) - i/„+i + 1] 



P 

p-i 



k=0 



p - k + 1 



)Sin,k) 



(3.4) 



p — k 



For a list of the resulting expressions up to p = 5 see section Q below. 

The following harmonic number identities can be expressed as a sum of this type: 



k=0 



fe=0 
P 



k=0 



For a list of the resulting expressions up to p = 5 see section Q below. 



4 Examples 



n 

fc=l 

n 

E k^"^-^ = Hj^^ 
fe=l 

" 1 



fc=0 



Formula ( 1.6 ): 



H,(r^' = ^n(n + l)(2n + l) 



(3.5) 



(4.1) 
(4.2) 
(4.3) 
(4.4) 

(4.5) 
(4.6) 

(4.7) 
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H^-'^^ - ^n{n + l){2n + l){3n^ + 3n - 1) 
H^-^'' = ^"2(n + l)2(2n2 + 2n - 1) 



Formula (2.14 1: 



J2 HkHa_k = {n + l)[(ff„+i - I f - hI^I, + 1] 

fc=0 
n 

kHkHr,-k = Hi^^[{Hn+i - If - H^^l^ + 1] 

fe=0 

k^HkH,,^k - - H^i] - ^n{n + l)(13n + 5)i/„+i 



fe=0 



+ ^"(" + l)(71« + 37) 

n 

^ PHkH^.k = - i^Si] - Y^n2(n + l)(7n + 5)H, 

+ ^n2(n + l)(35?i + 37) 



fe=o 



n+l 



1 



fe=0 



900 
1 



54000 

n 

1 



[n + l)(447n'^ + 468n^ + 17n - 32)i?„+i 
{n + 1) (20739^3 + 33066^2 + 4129n - 3934) 



fe=o 



360 
1 



21600 



n-'(7i + 1) (15771-^ + 218n^ + 17n - 32)i7„+i 

\n + l)(6839n3 + 1456671^ + 4129n - 3934) 



Formula (3.4): 



Y Hi = {n+ l)ff^+i - (2n + 3)i/„+i + 2(n + 1) 

fe=0 

n 

Y kHl = H^^^Hl^, - ir^+n- l)i/„+i +(„ + !)(„- 2) 



n 

Y^k'Hi ^ Hi-^^Hl+,--{4n'+9n^+5n+3)H^+i + — in+l){8n'+n+18) (4.19) 

A:=0 

n 

k'Hl = Ht'^Hl^, - -n(n + l)(n + 2)(3n + 
fc=o (4.20) 

l)(9n2 + 13n + 14) 

288 

n 
fe=0 

- — (72^5 + 31571"* + 410n3 + 13571^ - 32n - 30)i7„+i '^^'^"^^ 



900 



54000 

n 

5Z ^^^fe ^ ''^n+l 



(n + 1) (864n'* + 2241^^ + 262977^ + 1466?^ - 1800) 



fe=0 



—n(n + l)(7i + 2)(2n + 1)(10772 + 19n - 9)ff„+i '^"^'^^^ 
360 

^ (77 + 1)(20077'* + 73677^ + 115977^ + 971?7 - 366) 



21600 



Formula (3.5) 



J2 Hl-k = (" + ml+i ~ (2n + 3)iJ„+i + 2(77 + 1) (4.23) 

fc=0 

n ^ ^ 

kH^-k - Hl-^''Hl+i - ;^(3"' + 577 + l)i/„+i + -(77 + 1)(777 + 2) (4.24) 



k=0 



Y k^Hl-k = Hi-^'^Hl+i - ^(2?i + m^n^ + 17^ + 3)i?„+i 

fc=o (4.25) 

+ -^(77 + 1)(170772 + 10977 + 18) 

108 

n 

J2 k'Hl_, = Ht'^H^^, - -77(77 + l)(25r72 + 3777 + 10)i/„+i 

fc=o (4.26) 

+ + 1)(415772 + 42777 + 130) 

288 

71 

- ^(277 + l)(41l77* + 100277^ + 67977^ + 2877 - 30)i7„+i 

^ (n+l)(7211477'* + 10349177^+ 4612977^ + 146677 - 1800) 



54000 



fe=0 

- ■^n[n + l)(98n'' + 236^^ + 15971^ +n- 14)i7„+i '^^'^^^ 

+ n(n + l)(26978n'' + 49996^^ + 2959971^ + 196l7i - 3234) 

21600 

5 Computer Program 

The Mathematica® [5j program used to compute the expressions is given below, 
which should be added to the program listed in 

HarmSumR [0] = (n+1) ( (HarmonicNumber [n+1] -1) ~2-HarmonicNumber [n+1 , 2] +1) ; 
HarmSmnR [p_] : =HarmSumR [p] =Simplif y [1/ (p+1) 

(-n(n+l)~p HarmonicNumber [n+1] 

+Sum [Binomial [p , k] ( (n-k/ (p-k+1 ) ) HarmSumR [k] 

+(l+(-l)~k(p-k)/p(n+l)~(p-k-l))HarmSumF[k+l,l] 

+(p-k)/p HarmFun[-k-l]) ,{k,0,p-l}])] 
HarmonicSumR [p_] : =Module [{t=HarmTable [2] ,u}, 

u=Factor [CoefficientArrays [HarmSumR [p] ,t]] ; 

u[[l]]+Dot [u[[2]] ,t]+Dot [Dot [u[[3]] ,t] ,t]] 
HarmSumS [0] = (n+ 1 ) Harmoni cNumber [n+ 1] ~ 2- ( 2n+3 ) Harmoni cNumber [n+1] 

+2 (n+1) ; 

HarmSumS [p_] : =HarmSumS [p] =Simplif y [1/ (p+1) 

((n+1) ( (n+1) ~p(HarmonicNumber [n+1] '2-HarmonicNumber [n+1] ) 
-HarmonicNumber [n+1] +1) 

-Sum [Binomial [p , k] ( ( 1+k/ (p-k+1 ) ) HarmSumS [k] 

+ ( (p-k) / (k+ 1 ) - 1 ) HarmSumF [k+ 1 , 1 ] 

-(p-k)/p HarmFun[-k-l]) ,{k,0,p-l}])] 
HarmonicSumS [p_] : =Module [{t=HarmTable [1] ,u}, 

u=Factor [CoefficientArrays [HarmSumS [p] ,t]] ; 

u[[l]]+Dot [u[[2]] ,t]+Dot [Dot [u[[3]] ,t] ,t]] 
HarmSumT[p_] :=Simplify[ 

Sum[(-l)~k Binomial [p,k]n' (p-k) HarmSumS [k] ,{k,0,p}]] 
HarmonicSumT [p_] : =Module [{t=HarmTable [1] ,u}, 

u=Factor [Coef f icientArrays [HarmSumT [p] ,t]] ; 

u[[l]]+Dot [u[[2]] ,t]+Dot [Dot [u[[3]] ,t] ,t]] 



(* Compute some exaunples *) 
HarmonicSumR [3] //TraditionalForm 
HarmonicSumS [4] //TraditionalForm 
HarmonicSumT [2] //TraditionalForm 
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